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Abstract 

Extending a computation which appeared recently in pQ, we compute the transmis- 
sion and reflection coefficients for massless uncharged scalars and gravitational waves 
scattered by d > 4 Schwarzschild or d = 4 Reissner-Nordstr0m black holes, in the 
limit of large imaginary frequencies. The transmission coefficient has an interpreta- 
tion as the "greybody factor" which determines the spectrum of Hawking radiation. 
The result has an interesting structure and we speculate that it may admit a simple 
dual description; curiously, for Reissner-Nordstr0m the result suggests that this dual 
description should involve both the inner and outer horizons. We also discuss some 
numerical evidence in favor of the formulas of . 



1 Introduction 



Recently there has been a resurgence of interest in certain classical aspects of wave 
propagation in black hole backgrounds. Specifically, it has been suggested that scattering 
in the black hole background at large imaginary frequencies may carry information about 
quantization of the horizon area [SEISE], which has led to some recent investigations of the 
large-imaginary-frequency behavior. These recent studies [TIIdTTT] mostly focused on the large- 
imaginary-frequency "quasinormal modes," which are poles of the transmission coefficient 
and reflect the black hole's ringdown response to a perturbation, although 8\ also recently 
studied different aspects of the large-imaginary-frequency behavior using methods similar to 
those employed in PQ and in this paper. 

In Section |2] of this paper we present a computation of the transmission and reflection 
coefficients for waves traveling from spatial infinity to the black hole horizon, in the limit of 
large imaginary frequencies: uo — * +ioo or more precisely \ujth\ 3> 1, where r# is the horizon 
radius. We consider here the case of a massless scalar field; then for Schwarzschild black 
holes in d > 4, the result is 1 

P Pu _ i 

T(u)^— a , (1.1) 

R(u) w a 21 , (1.2) 
v > e^ + 3 K J 

and for Reissner-Nordstr0m in d = 4 we get 

e /3u> _ i 

T(uj) « — a =— , (1.3) 

v ; e^ + 2 + 3e"^ v ; 

R(u) ^iVs 5 5-. (1.4) 

Here j3 is the inverse Hawking temperature, and for Reissner-Nordstr0m (3i is the inverse 

u\r H 



Hawking temperature of the inner horizon; the symbol ~ means that there are O (l/ 
corrections both to the numerator and to the denominator. 2 The computation is a slight 
extension of one which recently appeared in [Q ; the new result is that we obtain the numer- 
ator as well as the denominator. We also calculate T(uj) and R(u) in the limit uj — ► —ioo. 
The results are summarized in Section l2~3l 

As we discuss in Section 12.41 the full result has a suggestive structure; it is strongly 
reminiscent of the results of piTU] , which showed that it is possible to reconstruct the effective 
SCFT description of near extremal black holes by studying T(u) for those black holes at 
small real frequencies. Essentially the idea is that the scattering of Hawking radiation back 
down the hole by the spacetime curvature modifies the blackbody spectrum in a precise way 

^^Thcre is a branch cut in T(u>) and R(u>) at positive imaginary frequencies; to switch to the other side of 
the cut one replaces uj by -uj in Ql.lj) . JQJ), (|1.4|) . 

2 In the case of Reissner-Nordstr0m the coefficient of these corrections becomes large in the Schwarzschild 
limit, which explains why the results do not simply reduce to those for Schwarzschild as /3/ — > 0. 
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and makes it agree with the spectrum obtained by an SCFT calculation. We are led to 
conjecture the existence of an effective description which determines the Hawking radiation 
at large imaginary frequencies and involves excitations with rather exotic statistics dictated 
by the denominators of (jl.3j) (such statistics were also proposed in for d = 4 

Schwarzschild). We discuss this conjecture in Section l2~4l 

Then switching gears, in Section |3] we discuss the numerical evidence supporting the 
analytical computations of [Tj for the denominators, some of which emerged only very recently 
[7j. The comparison to the numerical results also clarifies some subtleties about the relation 
between the formulas for Reissner-Nordstr0m and for Schwarzschild. 

We conclude with a few remarks on prospects for future progress. In two appendices, we 
discuss the analytic structure of the transmission and reflection coefficients and give a few 
more details of the asymptotic matching procedure used in the computation. 

2 The transmission and reflection coefficients 

In this section we will compute the transmission and reflection coefficients for waves 
propagating from infinity to the horizon of a static black hole. The analysis has been carried 
out so far for Schwarzschild in d > 4 and for Reissner-Nordstr0m in d = 4. To be concrete 
we mostly discuss Schwarzschild in d — 4; in Section l2~3l we give the results for other cases. 

2.1 The wave equation 

We consider wave propagation in the exterior region of a Schwarzschild black hole in 

d = 4, 

ds 2 = -f{r)dt 2 + /(r)" 1 ^ 2 + rW. (2.1) 
The coordinates ()2.1|) run from Th < r < oo, where 

r H = 2GM (2.2) 

and 

f{r) = l-r H /r. (2.3) 

The propagating field can be a minimally coupled massless scalar, an electromagnetic test- 
field, or a linearized perturbation of the metric; we label these three possibilities respectively 
by their spins j = 0,1,2. In each case the appropriate wave equation is separable. For 
example, in the j = case, writing 

<f>(r,t,0) =r^(r)Y lm (n)e iut (2.4) 

the massless Klein-Gordon equation A<p = becomes a one-dimensional wave equation on 
the interval — oo < x < oo, deformed by a potential: 

1>(x) = 0, (2.5) 



_ + V [r(x)} - cu 
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where 



V(r) = f(r) 



1(1 + 1) | r H 



(2.6) 



and we have introduced the "tortoise coordinate" x defined (up to an additive constant) by 

dr 



x(r) 



f(r) 



r + r H log(l - r/r H ). 



(2.7) 



Wave equations of the form (|2.5jl also characterize perturbations with j = 1,2 fUE]; only 
the form of V(r) changes. 3 The three cases can be nicely summarized, 



V(r) = f(r) 



(21 



In our computation it will be convenient to leave j generic and then take the limit j — > or 
j — > 2 in the answer. We note that our analysis will depend on the 1/r 4 term in V(r) being 
dominant; it is therefore singular as j ' ^ 1 and we cannot predict the behavior there. 



2.2 Transmission and reflection in d = 4 Schwarzschild 

In this section we will compute the transmission and reflection coefficients for ([2.5)1 . in 
the limit of large imaginary frequencies. Most of the ideas required to do this calculation 
have appeared already in P , in particular the idea of using the boundary condition at Th to 
provide the monodromy. The general method of analytic continuation to complex r is much 
older and has been used successfully for numerical calculations in the past (see [H] for a 
recent discussion and review.) 

Although we are studying a wave equation defined on the physical region r# < r < oo, to 
get the analytic continuation to imaginary u we have to consider complex r. The behavior 
of the tortoise coordinate ()2.7|) at complex r plays an important role in the computation; in 
particular, since x(r) contains r?j log(l — r/rs), it is multivalued, with the branches differing 
by Ax = 2irinrH- This multivaluedness is crucial for our computation. Re x, however, is 
uniquely defined; in Figured] we have indicated its sign by shading. 

The idea of the computation is to study the behavior of a solution ip(r) to (|2.5|) as we 
travel around a specific contour 7, shown in Figure ^ We will calculate the effect of the 
trip around 7 in two different ways: a "local" computation by integrating the differential 
equation directly and a "global" computation using the boundary conditions at the horizon. 
Comparing the two gives enough information to calculate the transmission and reflection 
coefficients. 

Since V(r) vanishes faster than l/\r\ as r — > 00, as we run to infinity in any direction 
the asymptotics of a solution to (|2.5J) are of the form 

V>(r) ~ C + e +iujx + CLe"^. (2.9) 

3 For j = 2 there is an added complication because one can consider either "axial" or "polar" perturbations; 
(12.81) strictly applies only to the axial case, but the two are closely related and have the same quasinormal 
frequencies |13| . 
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Figure 1: The complex r-plane. The Re (x) < region is colored gray, and the points r = 0, 
r = th and a contour 7 for calculation of the transmission and reflection coefficients are 
marked. 7 runs along the line Re (x) = to \ux\ ^> f in each direction. The symbol + 
represents the direction in which the boundary conditions "at infinity" are imposed. 

In (12. 9 j) there is an ambiguity in defining e ±lu;x , which arises from the multivaluedness of 
x(r). This ambiguity implies a compensating ambiguity in defining C±, even for a fixed 
ijj(r). This ambiguity is compounded by the fact that ip(r) itself can have multiple branches, 
since it is the solution to a differential equation which has singularities. The upshot is that 
as we travel around 7 and return to our starting point the coefficients C± do not return 
to their original values (i.e. they have a nontrivial monodromy), because we end up on a 
different branch both for x(r) and for ip{r). A 

As explained in Appendix [21 there is a branch cut in T(u) and R(u) on the positive 
imaginary tu-axis. If we are on the right side of the branch cut (taking Re u infinitesimally 
positive) then T(u) and R(u) are defined by the boundary conditions 



where ip{r) is a solution of (|2.5|) . As was done in PQ, we will obtain our results by computing 
in two different ways the monodromy of the coefficient of e~ tuJX . 

To make this last statement clearer, let us study for a moment the behavior of ip(r) if we 
simplify by setting V(r) = 0. In this case we may replace the asymptotic relation in (|2.10j) 
by exact equality: ip(r) = e tuJX + R(uj)e~ luJX near -j*. Then we follow ip(r) along 7. Since 

4 Actually, our boundary conditions at rjj will imply that only C_ has nontrivial monodromy. 




as r 



as r 



* 



(2.10) 



locally we always have free wave propagation in the variable x, when we return to ^ we still 
have ip(r) = e tuJX + R{uo)e~ lwx , i.e. the coefficient of e~ luJX is the same before and after the 
trip around 7. On the other hand, 7 encloses a singularity r = around which both x(r) 
and ip{r) are multivalued. The boundary condition (j2.10|) implies that ip(r) picks up a phase 
e 2nujr H on an y pg^k encircling rn in the clockwise direction. Furthermore, e~ luJX is multiplied 
by the phase e~ 2nujTH on such a path. These two phases together imply that the coefficient 
of e~ tuJX must be multiplied by e 4nujra after the trip around 7. But we already showed that 
it is unchanged after this trip; so we get R(u) = e 47IWj? R(co) , which implies R(u) = 0. 
This is as expected: there is no reflection, since there is no potential! We can calculate 
T(u) by traveling on any contour running from ^ to the horizon; then, since we have free 
wave propagation, the boundary conditions at it given by (j2.10J) are propagated directly 
to the horizon, and comparing with the boundary condition there we see immediately that 
T(u) = 1. (We also see immediately that R(u) = 0, which we had already concluded above 
by more elaborate means; the reason for the more complicated discussion above is that it 
generalizes to the case of a nontrivial potential.) 

Now let us consider the real problem at hand, where V(r) is given by ()2.8j) : 

K (r ) = /(r)p^ + ^). ,,U) 

The arguments above still show that R is multiplied by the monodromy e Anujr ' H after a trip 
around 7; this depended only on the boundary condition (j2.10j) at r# and the form of x(r). 
Note that in the black hole context Attth is the inverse Hawking temperature j3 = 1/T#, so 
the monodromy can be written e^. 

Now we follow 7 beginning at where we have the boundary condition (J2.1(J|) . In 
Appendix A we describe a natural basis ip±(r) of solutions near r = 0; as discussed there, by 
asymptotic matching along the line Re x = we can relate the form of if) slightly southwest 
of r = to the form of ip & t Namely, if we write the solution southwest of r = as 
if)(r) = A + ip + (r) + A-i[)-.(r), and introduce the notation 

(m) = A + e im{1+3)w/ * + A_e m{1 - j)7T/ \ (2.12) 

then matching the asymptotics (jB.10|) of tp± to ()2.10j) gives 

(-l)wl, (2.13) 
(l)wfl(u;). (2.14) 



Here ~ means there are O (l/ \J\lu\th] corrections. Next we want to express T(lj) in terms 
of A±. For this we first make a quarter-turn from southwest to southeast, then travel out 
along Re x = until \ux\ ^> 1, and then run directly to r = Th- Since on this trip Re x < 
and the u 2 term dominates V(r), the WKB (physical optics) approximation guarantees that 
the coefficient of e iuJX will not be modified as we travel to the horizon, to leading order in uj. 5 

5 We are being slightly loose here, but the idea is that e luJX dominates the ui — > ioo asymptotics when 
Re x < 0, so any modification to its coefficient would cause those asymptotics to deviate from the physical 
optics approximation. This is the same idea used in the old derivations of the WKB transmission formula 
for a particle with energy near the top of a barrier [15j : we will use it several times. 
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Therefore this coefficient at the horizon is determined simply by the asymptotics of tp as we 
run southeast from r = 0, which are given by (|B.8|) and (jB.10|) . leading to 



(3)«T(w). (2.15) 

Finally we want to complete our trip around the contour. So we make a three-quarter 
turn from southwest to northwest, then travel north until \ujx\ ^> 1. Then the asymptotic 
matching using (jB.8|) . (jB.lO)) determines the coefficient of e~ tuJX to be (5). On the other 
hand, we can follow this coefficient around the big semicircle and back to again by the 
WKB approximation it must be unchanged by this trip. But as discussed above, after the 
full round trip around 7 the coefficient of e~ luJX is e /3w i?(u;), so we get 



R{uj). (2.16) 



The relations (J2.13|) . (|2.14j) . (J2.15|) . (|2.16|) define an inhomogeneous linear system for the 
four unknowns A + , A_,T(uj), R(uj). Solving it gives the final answer: 



T(w) ~ —a 7 r, (2-17) 

v ; e#" + (l + 2cos7rj)' v ; 

, v COS ^7 

R(uj) « 2i— -. — ^ -. 2.18 

v ; + (1 + 2cos7tj) v ; 



Here the symbol ~ means there can be O y-/ \j \^\ r Hj corrections both in the numerator 
and the denominator. The previous analytical results of PjJIH] for the asymptotic quasinor- 
mal frequencies are recovered by looking at the poles of R(uj) and T(uj), which occur at uj 
satisfying 

+ (1 + 2C0S7T?) « 0. (2.19) 

In particular, for j = or 2 we recover Re to ~ Tjj log 3, which was important for 00] • 

The results (|2.17|) . ()2.18|) apply for u> approaching the branch cut from the right, i.e. Re uj 
infinitesimally positive. If instead we took Re uj infinitesimally negative then the boundary 
condition at infinity in ()2.10|) would be set at the north end of Figure Q rather than at -fa, and 
7 would run in the opposite direction. The effect is to replace e^ w by in (|2.1fi|) and hence 
to replace uj by — uj in the final result (still valid for uj — > zoo.) Then we recover for j = or 
2 the quasinormal frequencies on the other side of the imaginary axis: Re uj ~ —Th log 3. 

We can also study T(— uj) and R(—u) as uj — > zoo. The analysis is similar to that given 
above with two differences. First, there are some signs that have to be reversed in the 
asymptotics (jB.lOJ) . Taking these signs into account, the analogues of (j2.13j) . (j2.14)l . (|2.15jl 
are 

(1) w 1, (2.20) 

(2.21) 

(l)wT(^w). (2.22) 
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Second, the WKB approximation on the large semicircle is now good for the coefficient of 
e %uix j-g^hgj. than for e~ luJx . For this coefficient the two monodromy factors e 27TWTH cancel 
instead of multiplying, giving instead of (|2.16j) 

(5) sa 1. (2.23) 

Solving, we obtain 

T(-co) « 1, (2.24) 

R(-uj) « -2i cos -irj. (2.25) 

One can also see that T(— uj) ~ 1 directly just by noting that (as discussed in Appendix lA"|) 
when Im uj > the analytic continuation is not necessary for defining T{—uo) and R(—uj); we 
can just work directly on the real line, and then the WKB approximation gives T(— uj) m 1 
immediately since there are no turning points. 

2.3 Transmission and reflection for other black holes 

The computation of the transmission and reflection coefficients given in Section EI21 gen- 
eralizes in a straightforward way to the other black holes considered in jTj, with modifications 
in each case just like the modifications made in PQ for the computation of the quasinormal 
frequencies. Here we only report the answers. For Schwarzschild in d > 4 coupled to a 
massless scalar 6 the result is just as it was in d — 4: as uj —>■ zoo, 

e /3w _ 1 

T(u) « — a , (2.26) 

v ; e^ + 3' V ; 

R(uj) « — — , (2.27) 
v i e^ + 3 y ' 

T(-uj) « 1, (2.28) 
R(-u) « -2i. (2.29) 

For Reissner-Nordstr0m in d — 4 coupled to massless scalar or electromagnetic-gravitational 
waves the result is 

T(u) « -3— = — s — , (2.30) 

v ; e^ + 2 + 3e"^ v ; 

1 + e~^ IU 

R(u) « izv^-^ 1 3—, (2.31) 

r(-w) « 1, (2.32) 
« Tiv^, (2.33) 



6 The restriction to scalar fields arises because the appropriate generalization of (|2.8|) is not known for 
arbitrary metric perturbations in d > 4. 
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where the sign is determined by the type of field we consider. A useful check on the results 
is provided by the analytically-continued conservation of flux, 



T(u)T(-u) + R(u)R(-u) = 1. 



(2.34) 



2.4 Hawking radiation and greybody factors 

Taken together, the results of Section 12.31 for the transmission coefficient at large imagi- 
nary frequencies have a rather suggestive structure. To understand it, let us first recall one 
physical role played by the transmission coefficient at real frequencies. 

The spectrum of Hawking radiation seen by a static observer at spatial infinity is not 
thermal; rather, the particle flux obeys [To] 



The interpretation of this formula is that Hawking radiation is produced with a thermal 
spectrum at the horizon and then the spacetime curvature between the horizon and infinity 
scatters some of the radiation back down the hole, so that the thermal spectrum is multiplied 
by the "greybody factor" |T(u;)| 2 . 

In a remarkable paper |Hj it was shown that, for a certain class of near extremal five- 
dimensional black holes admitting a D-brane construction jT2j, and for energies small com- 
pared to 1/th, the greybody factor allows semiclassical general relativity to reproduce the 
emission spectrum of the appropriate system of D-branes. Namely, for small u and in a 
certain "dilute gas" regime of the black hole parameter space, a direct classical computation 
gives 



Here (3l and (3r are parameters in the classical black hole solution, which can be identified 
with inverse temperatures characterizing the left-moving and right-moving excitations in an 
effective SCFT which describes the D-brane degrees of freedom. The factor e® w — 1 in ()2.36|) 
cancels the denominator of ()2.35|) . leaving the spectrum observed from infinity: 



But ()2.37j) is exactly the emission spectrum calculated from the effective SCFT! To quote 
"to the observer at infinity the black hole, masquerading in its greybody cloak, looks like the 
string, for energies small compared to the inverse Schwarzschild radius of the black hole." Of 
course, this result was discovered after the D-brane construction was already known; but if 
the order had been reversed, one might have speculated that this emission spectrum gave a 
potential clue to the microscopic description of the black hole. This perspective was pursued 
further in ^0] where the calculation of the greybody factor was extended to near extremal 
charged rotating black holes, both in d = 4 and d = 5, and again found to be consistent with 
an effective SCFT description. 




(2.35) 




(2.36) 




(2.37) 
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In the present case of the Schwarzschild or Reissner-Nordstr0m black hole we have a 
similar result, with the important caveat that it does not apply directly to the emission 
spectrum at real uj. Rather, suppose an observer at spatial infinity measures the exact 
emission spectrum (for a massless scalar field minimally coupled to the curved background) 
over some interval of real uj and then analytically continues to large imaginary uj. What is 
the result? Using the fact that T{u)* = T(-u) for real uj, together with (jQHj) and (l2~2"8l . 
the analytic continuation of |T(a;)| 2 to large imaginary uj is 

T(w)T(—u) « (2.38) 

So substituting in ()2.35|) we find the same cancellation between the numerator and denom- 
inator that occurred in P; the observer at infinity would compute a flux which, for large 
imaginary uj, depends on uj as 

F(uj) oc . (2.39) 

Similarly, in the case of the Reissner-Nordstr0m black hole the result is 

F(u) oc -5 -— , (2.40) 

where /3j is the inner horizon temperature. 

What could these formulas mean? The fact that they are always written in terms of Boltz- 
mann factors (or more properly "Boltzmann phases" since uj is pure imaginary) is suggestive 
and we propose that, just as the spectrum of Hawking radiation of a near extremal black hole 
at small real uj can be obtained from an effective SCFT, the spectrum of a Schwarzschild or 
Reissner-Nordstr0m black hole at large imaginary uj also admits an interpretation in terms 
of some new degrees of freedom. These degrees of freedom would be expected to have rather 
exotic statistics (the possibility that the quasinormal modes of d = 4 Schwarzschild could 
imply "tripled Pauli statistics" appeared also in 6 .) In the Reissner-Nordstr0m case the 
occurrence of both e@ w and e@ lUJ suggests further that they should live on both the inner and 
outer horizons (such descriptions have been considered before, e.g. in [TH].) At any rate, 
the bulk scattering amplitudes we have computed could then be related to correlation func- 
tions in this conjectural "boundary" theory. Since our result is valid around large imaginary 
frequencies it is natural to suspect that this boundary theory bears some relation to the 
Euclidean continuation of the black hole. 

One might wonder why the structure we observe was not already visible in ^Uj. After 
all, that paper computed precisely the same function we are computing. The reason our 
T(uj)T{— uj) is not simply the analytic continuation of the |T(a>)| 2 in ^U] is that both are 
approximations to the true answer, valid in different regimes. If we know a function f(oj) 
exactly at real uj, then we can analytically continue it to get the exact answer for imaginary uj; 
but if we neglect a small correction at real uj it can invalidate the continuation to imaginary 
uj. This problem also occurs when using AdS/CFT to obtain Minkowski space boundary 
correlators, as discussed recently in [T^ OUj . 
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3 Comparison to numerical results 



Since some of the asymptotic formulas computed here and in (HE] are still new and 
the result for Reissner-Nordstr0m is rather strange, it may be useful to summarize the 
numerical evidence supporting them. While I am not aware of any numerical results on the 
full asymptotic transmission coefficient, there have been several studies of the asymptotic 
quasinormal frequencies for various types of black hole, to which the analytical results may 
usefully be compared. 

In the case of the Schwarzschild black hole the only numerical results are in d — 4 and 
they support the conclusion that the asymptotic real part of the frequency is T^log3, at 
least for j = 2 perturbations |2*T1|22|I231I24| . I do not know of any literature on the asymptotic 
quasinormal frequencies for other j. In (22! the corrections at finite uj were also tabulated; 
the results there are consistent with the hypothesis that the first subleading correction is 



O fl/y Mre") an d proportional to (I 2 + I — l). 7 This is consistent with an estimate of 

the correction to the Bessel function asymptotics as discussed in Appendix |Bl but I have 
not managed to evaluate the exact coefficient to compare with [22] . Recently jH] appeared 
containing a computation closely related to the one given here, obtaining the discontinuity 
across the branch cut of the solution satisfying the outgoing boundary condition at infinity 
as well as the first correction; this correction was found to be proportional to (/ 2 + / — 1) 
and it is likely that those results can be used to get the correction to the transmission and 
reflection coefficients as well. 

In the case of the Reissner-Nordstr0m black hole in d — 4, the asymptotic formula for 
the quasinormal frequencies is more intricate. It is convenient to introduce the parameter 



where rj and r# are the inner and outer horizon radii respectively; a runs from 1 for the 
uncharged black hole (Schwarzschild limit) to oo for the extremal black hole. In terms of a 
and uj = 8nGMu, the requirement that the transmission coefficient (jl.3j) has a pole becomes 



The behavior in the asymptotic limit depends crucially on whether a is rational. If a = p/q 
then we can write x = e^ 9 and (|H.2|) becomes a polynomial equation in x, with finitely 
many solutions. After taking the logarithm each solution then gives an infinite family of 
quasinormal frequencies, with constant real part and evenly spaced imaginary part; so instead 
of just one asymptotic real part we have a finite number of them. For example, if a = 2 we can 
set x = e u and then get x 2 + 3x + 2 = 0, which has two solutions x = —1, x = —2, giving two 
families of asymptotic quasinormal frequencies, at uj — 2ni{n + ^) and uj = 27ii(n + 1) + log 2. 
On the other hand, if a is irrational, the behavior of the asymptotic quasinormal frequencies 
is more complicated; it probably admits some statistical description in terms of the continued 
fraction representation of a. Such a sensitive dependence on whether a is rational seems 

7 This proportionality was discovered by Martin Schnabl. 





+ 2^0. 
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Figure 2: Inspiraling behavior of the asymptotic formula (|3.2|) for the n = 5 quasinormal fre- 
quency as the charge is increased toward extremality. The axes are Re GMu and Im GMu. 
The center of the spiral is at GMu = (ll/8)i + log3/87r. 

strange from the point of view of general relativity, but it might not be so strange from the 
point of view of an underlying microscopic description of the black hole which makes the 
parameters discrete. 

We can also follow the behavior of a single quasinormal frequency as we increase a from 
1 to oo. At a = 1 the solutions of (|3.2|) are simply u = 2ni{n + |) + log 5. 8 Increasing a 
continuously the frequency traces a spiral in the complex plane, shown in Figure 121 for n = 5. 
The center of the spiral corresponds to the extremal limit a — > oo in which we may neglect 
the last term of (|3.2j) . obtaining simply e 1 ^ + 3 ~ 0, which implies 



Amusingly, this result in the extremal limit agrees with the asymptotic quasinormal fre- 
quency for a Schwarzschild black hole of the same mass. 

The spiral structure shown in Figure El was visible already in the low-order modes com- 
puted in [23]. But very recently the full asymptotic formula has been confirmed in detail 
by the numerical analysis of [7], and looking at their results is useful not only to convince 

8 The reader may be alarmed that a = 1 is the Schwarzschild limit and we are getting Re u> w Tjj log 5 
instead of the correct Re w « Tjj log 3 for Schwarzschild. Such a reader is urged to suspend disbelief for a 
few more paragraphs. 




(3.3) 
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oneself that the formula is asymptotically correct but also to understand the nature of the 
0(l/ \J\u\rnj corrections. Rather than looking at the spiral it is convenient to plot Re uj and 
Im uj separately as functions of the charge Q; this is done in Figure|3]for n = 30, 5000, 10000. 9 
We see that at large n the numerical and analytical results agree extremely well, except 
at small Q where the real part is significantly different even at n = 10000. How is such 
a disagreement possible in the asymptotic limit? The point is that, as already mentioned 
in pp, for the asymptotic formula to be valid one has to take the n — * 00 limit at fixed 
Q 7^ 0. In this limit we will eventually find agreement with the asymptotic formula (e.g. if 
we fix Q = 0.16, we can see from Figure El that n = 5000 is not sufficiently large for the 
asymptotic behavior to take over, but n = 10000 is.) The required n diverges as Q — > 0, 

or put another way, the O (l/ ^J\oj\rHj corrections to the asymptotic formula (|3.2j) appear 
with a coefficient which diverges as Q — > 0. A more systematic analysis of these corrections 
should demonstrate this behavior explicitly; here we only remark that the analysis in pQ 
which led to ()3.2|) is certainly singular in the limit Q — > 0, so we would expect that the 
corrections would indeed blow up. 



4 Conclusions 

In this paper we have presented some results about the scattering of free massless quan- 
tum fields by a black hole, in the limit of large imaginary frequencies. These results stand 
independent of any speculation about their meaning, but their form is tantalizing and it 
seems that they deserve to have a simple explanation. The remarks about a possible bound- 
ary interpretation which appear in Section E31 are of course pure conjecture at the moment; 
hopefully the situation will be clarified in the future. 

It remains unclear whether one can connect the computation presented here to the area- 
quantization ideas which motivated the consideration of the high-imaginary-frequency quasi- 
normal modes in j2H2|. It is clear that these ideas will have to be modified or reinterpreted 
in some way if they are to be correct, since the Reissner-Nordstr0m and Kerr black holes 
do not share the simple behavior seen in the Schwarzschild case. Nevertheless, it is possible 
that this could be done, and the spin-network picture discussed in |3] could be a candidate 
for the boundary description we are conjecturing. 

There are several other calculations along the same lines which could be done and might 
shed some light on the situation. One possibility would be to include chemical potentials, 
either by studying a rotating black hole or a charged scalar in the Reissner-Nordstr0m 
background. In both cases one faces the technical difficulty that the potentials appear to be 
long-range (they involve 1/r) and hence the behavior is not simply plane- wave at infinity. In 
the rotating case, at least for the Kerr black hole, this problem can be avoided by a clever 
change of variable introduced by Sasaki and Nakamura [2H], 10 but the problem remains 

9 I am indebted to E. Berti and K. Kokkotas for sharing their data files with me, and to L. Motl for 
preparing earlier versions of these figures as well as suggesting that they be made in the first place. 
10 I thank Shinji Mukohyama and Eanna Flanagan for informing me about this reference. 
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Figure 3: Real parts (left) and imaginary parts (right) of quasinormal frequencies for the 
Reissner-Nordstr0m black hole, at mode number (top to bottom) n = 30, 5000, 10000. The 
red line is the numerical data of [7j and the black line is the asymptotic formula (|3.2|) . The 
vertical axis is (Re or Im) GMuj; the horizontal is Q, in units where Q — 1 is extremal. 
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difficult because the angular equation also has to be treated approximately (the centrifugal 
potential, which was just /(/ + l)/r 2 for Schwarzschild and had no effect to leading order in 
u, depends on u for Kerr.) For Reissner-Nordstr0m I am not aware of a change of variables 
which makes the potential short-range, but it should be possible to make one and this 
problem might be more tractable than Kerr. On the other hand, in the Kerr case there are 
some numerical results available [7j which look puzzling but might be helpful. One could also 
study the high-imaginary-frequency behavior for fermionic fields. Here the relevant radial 
equation has been given e.g. in [27], and the low overtone quasinormal modes were recently 
considered in |2~Hj . 

Another possibility would be to apply our method to the five-dimensional near extremal 
black holes studied in Plll7j . for which we already have not only an effective SCFT but a full 
microscopic description, and see whether the resulting large-imaginary-frequency behavior 
can fit into that description somewhere. 
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A Generalities on one-dimensional wave equations 

The notions of transmission and reflection coefficient for a one-dimensional wave equation 
are well known, at least for real frequencies. However, there are some subtleties associated 
with analytic continuation to complex frequencies. Under fairly general conditions this 
analytic continuation does exist; the questions are what its singularities are and whether we 
can construct it in a way which facilitates explicit computation. 

Suppose given a wave equation on the real line — oo < x < oo with the form of (|2.5|) . 



where V(x) — > faster than l/\x\ as x — > ±00. We can then define T(a>), R(u) for waves of 
real frequency uj 7^ by fixing the boundary conditions 



dx 2 



+ V(x) - co 2 4){x) = 0, 



(A.l) 




(A.2) 
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The boundary conditions (jA.2j) also suffice to extend T(u) and R(u) to complex-analytic 
functions on the lower half-plane Im uo < 0. Furthermore these functions are nonsingular; 
there cannot be any poles because they would imply the existence of a normalizable solution 
to (jA.l|) with complex uo, which is impossible because of the self-adjointness of —d 2 /dx 2 + 



It is not straightforward to define T{uj) or R(uj) for uo in the upper half-plane, because 
in that case the boundary condition (jA.2J) as x —>■ +00 does not give a constraint on the 
solution; the solution e +tuJX is exponentially small compared to e~ tuJX as x —>■ +00, so adding 
e +iux with an arbitrary coefficient does not affect the asymptotics asx-> +00, and there is 
no way to define what we mean by requiring this coefficient to be 1. In some cases we can 
nevertheless analytically continue T{uo) to the upper half-plane; the existence and singularity 
structure of such a continuation depend on the properties of V(x). We now turn to some 
examples. 

The simplest possibility is V(x) = for all x. In this case the prescription above gives 
simply T(uj) = 1, R(uo) = for all uo in the lower half-plane, which has the obvious analytic 
continuation T(uo) = 1, R(uo) = for all uo. 

Another possibility is that V(x) 7^ but V(x) strictly vanishes as x — > ±00 (e.g. a 
finite square well). In that case, for any uj, the solutions of (jA.l|) near ±00 are strict linear 
combinations of e ±luJx . Therefore T(uj) can be uniquely defined for any uj just by replacing 
the asymptotic relation ~ by exact equality in (|A.2|) . This prescription defines T(uo) and 
R(uo) univalently for any u> 7^ 0, so T{uj) and R(uo) are single-valued, although they can have 
isolated singularities at uo = and in the upper half-plane. 

A third possibility is that V(x) admits analytic continuation to complex x and further- 
more V(x) — > faster than l/\x\ as x — > 00 in any direction (e.g. V(x) = l/(x 2 + 1).) Note 
that this case and the one just considered are mutually exclusive. In this case we can define 
T(uo) by setting our boundary conditions on the line uox G R rather than 



The advantage of (jA.3|) over (jA.2|) is that the solutions are oscillatory on the line uox G R and 
hence the asymptotics (|A.3|) determine T(uo) and R(uo). There is one additional complication: 
if V(x) has singularities along the line uox G M, then i()(x) will generally be multivalued, so 
we have to specify the contour along which we travel from one side of the line uox G R to 
the other. To make T(uo) and R(uo) holomorphic in uo the only option is to continuously 
deform the contour by "Wick rotation" from the real axis; but then after a full 2tc rotation 
the contour will hang on the singularities, so T(uo) and R(uo) are potentially multivalued. 

Finally we come to the case which actually occurs for black hole transmission coefficients, 
say for d = 4 Schwarzschild. Here the analytic structure of T(uo) and R(uo) has been exten- 
sively studied because it is relevant for the long time behavior of black hole perturbations; 

n This self-adjointncss fails for non-normalizable ip(x), which formally explains the existence of quasinormal 
modes, i.e. poles of T(uj) and R(uj) in the upper half-plane. 



V(x). 
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(A.3) 
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see in particular In this case V(x) = V[r(x)], with r(x) defined implicitly by (|2.7|) . 

and while V(r) is manifestly single- valued all over the complex r-plane, r(x) generally has 
multiple branches (because of the logarithm in (|2.7]l ). So the analytic continuation of V(x) is 
a multivalued function of x. Nevertheless, since the only singularity of V(r) occurs at r = 0, 
the singularities of V[r(x)] can only occur at x for which some branch of r(x) has r(x) = 0, 
namely at x — 2mn. So T(u), R(uj) can be defined by Wick rotation of the x-plane contour 
to any uj with Re to 7^ 0. However, at Re uj = we encounter an infinite chain of singulari- 
ties and we have to choose which way to go around them; our choice will depend on which 
way we made the Wick rotation (i.e. whether we are coming from Re uj > or Re uj < 0). 
T(u), R(u) therefore have branch cuts at positive imaginary uj. In the body of the paper we 
compute these functions for uj on the right side of the cut, i.e. Re uj infinitesimally positive; 
in that case, after performing the Wick rotation in the x-plane, (|A.3|) corresponds to (see 
Figure □ in Section l2~2J) 

|7V)e~« (A.4) 
\e luJX + R(uj)e~ luJX asr^*. 

If we had assumed Re u < instead, we would have put at the north of Figure [T] instead 
of the south. 



B Asymptotics of solutions to the wave equation 

In our analysis of the solutions of (|2.5|) along the line Re x = we need to know their 
asymptotics in the limit uj — > zoo. In this limit the term uj 2 dominates the potential every- 
where except near r = 0, so we expect that the interesting behavior will occur there. Near 
r = we have from (|2.7J1 

x /r H = r/r H + log(l - r/r H ) = ~^(r/r H ) 2 - ^(r/r H ) 3 + O ((r/r H ) 4 ) (B.l) 

which can be inverted to obtain 

r/r H = ±i^2x/r H + 2x/3r H + O ([x/r H f /2 ) (B.2) 

Substituting this in (|2.5jh (|2.8J1 we obtain the leading-order behavior of ()2.5|) near the r = 
singularity together with its first correction: 



d 2 j 
+ - 



dx 



Ax 2 



±i- 



V 



x 



3/2 



+ 0((xr H ) 



UJ 



V = 0, 



with 



3/(/ + l) 



6V2 



(B.3) 



(B.4) 
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Further setting p = oox, 



dp' 



+ 



±i- 



V 



p 3 / 2 y/ur H 



+ 



ip = 



(B.5) 



Integrating the subleading terms twice in p shows they can at worst give corrections where 
Si/j/ip scales as a power of p/uoru- Then in the u — > ioo limit with p 1, one may neglect 
this term as well as the constant term 1 and write simply 



d 2 



+ 



if) = 0. 



dp 2 4p 2 

So when p 1 there are two possibilities for the asymptotics: 



tp(p) ~ p 



(i±j)/2 



(B.6) 



(B.7) 



It turns out that both possibilities occur and we can pick two linearly independent solutions 
ip± by fixing their rotation properties (at least for non-integer j): 



On the other hand, for 1 p, (jB 



is simply 

d 2 



dp' 



+ 1 



■tp = 



(B.8) 



(B.9) 



and so the solutions are asymptotic to free waves e ±%p . The issue is now to interpolate 
between these two regions. To leading order in u this can be done by neglecting the \j \/To 
term in (jB.5|) ; one then finds that ip± are asymptotic to y/pJ±j/2(p) in the region p <^ \w\rn, 
and since this region overlaps with the region 1 <C p in the limit of large u, one can match the 
known asymptotics of Bessel functions to the free wave propagation. This was the approach 
taken in pQ and it gives for 1 p 



MP) ~ cos(p - (1 ± j>/4) + O l/J\uj\r H 



;b.io) 



From ()B.4|) . ()B.5|) we see that the O (l/ ^J\u\r H ^J correction to ()B.10|) should be proportional 
to (j 2 - 1) - 3/(7 + 1); in the case j = 2 this dependence was also found in where in 
addition the exact numerical coefficient was evaluated. 
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